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Abstract. We report BRAHMS results from RHIC d+Au and p+p collisions 
at y/SNN = 200GeV. A remarkable change in the nuclear modification factor 
RdAu is seen as the pseudorapidity of the detected charged hadrons changes 
from zero at mid-rapidity to 3.2 at the most forward angle studied during 
the 2003 run. For pseudorapidity rj > 1 the suppression of the R cp factor is 
more pronounced in the sample of central events in contrast to the behavior at 
mid-rapidity where the central events show higher enhancement compared to 
a semi-central sample. These results are consistent with a saturated Au wave 
function strongly affected by quantum evolution at higher values of rapidity. 
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1. Introduction 



j_i ■ The first BRAHMS results from d+Au collisions at forward rapidities have gen- 

erated heated discussions since they were presented in their preliminary form at 
the DNP meeting in Tucson AZ. At that time, the theoretical work offered two 
clearly differentiated views. On one side, the groups that postulate the formation 
of the Color Glass Condensate CGC at RHIC, had results that demonstrated 
the presence of Cronin enhancements [|2j in CGC 0J |SJ as well as studies 
that included quantum evolution to describe how the nuclear modification factor 
would be modified as the collisions are studied at higher rapidities [ [7J |Hj and [ 
9 . These groups describe the d+Au collisions as coherent multiple interactions 
between the deuteron valence quarks and a saturated Au wave function at small 
values of x. As the density of gluons in the Au nuclei grows, higher order correc- 
tions to the gluon density are included within the formalism referred as quantum 
evolution. The net effect of these corrections is an overall reduction of the number 
of gluons compared to scaled p+p collisions. The effect of quantum evolution is also 



2 



R. Debbe et al. 



present in the centrality dependence; the more central the collisions, the stronger 
the effect of quantum evolution making the suppression of the number of gluons 
more pronounced in central events. 

Other groups had worked the problem based on an standard description of 
the Cronin enhancement as incoherent multiple scattering at the partonic level [ 
1101 1111 I12j . This description of the d+Au collisions would continue exhibiting 
a Cronin type enhancement close to the deuteron fragmentation region, and the 
strength of the enhancement would increase with the centrality of the collision. 

We present here a description of the BRAHMS results [EH within the context 
of saturation that includes the effects of quantum evolution. We are well aware that 
even though our results are consistent with the descriptions offered by the theory 
that includes the presence of saturation in the initial state, more measurements are 
necessary to eliminate other explanations. 

2. Experimental results 

The spectra presented in this contribution were extracted from data collected with 
both BRAHMS spectrometers, the mid-rapidity spectrometer (MRS) and the front 
section of the forward spectrometer (FFS). A detailed description of the BRAHMS 
experimental setup can be found in [ I14| . The low multiplicity of charged particles in 
the proton+proton and d+ Au collisions required an extension of the basic apparatus 
with a set of scintillator counters (called INEL detectors) . These detectors [ HK| 
cover pseudo-rapidities in the range: 3.1 <| r] |< 5.29, and define a minimum biased 
trigger. This trigger is estimated to select w 91% ± 3% of the 2.4 barns d+Au 
inelastic cross section and 71% ±5% of the total inelastic proton-proton cross section 
of 41 mb. The INEL detector was also used to select events with collision vertex 
within ±15 cm of the nominal collision point with a resolution of 5 cm. 

The centrality of the collision was extracted from the multiplicity of the event 
measured within the angular region | 77 |< 2.2 with a combination of silicon and 
scintillator counters [ llfij . 

2.1. Spectra 

Figure^shows the invariant yields obtained from p+p collisions (panel a) and d+Au 
collisions (panel b). For each system we studied particle production at 40 degrees 
with the MRS spectrometer and 12 and 4 degrees with the FFS spectrometer. Each 
distribution was obtained from several magnetic field settings and corrected for 
the spectrometer acceptance, tracking and trigger efficiency. No corrections were 
applied to the spectra for absorption or weak decays. Statistical errors are shown 
as vertical lines, and an overall systematic error of 15% is assigned to each point. 
The p+p spectra have also been corrected for trigger efficiency by 13 ± 5% to make 
them minimum biased with respect to the total inelastic cross section. We fitted 
the spectra at 77 = 3.2 with a power law function . t+ . n and the integral of that 
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function over pj, is compared for consistency in table ^ with UA5 results [ El f° r 
the p+p system, as well as our own multiplicity measurement for the d+Au system 

[El- 




Fig. 1. Spectra for charged hadrons at different pseudo-rapidities. Panel a shows 
the spectra obtained from proton-proton collisions and panel b those from d+Au 
collisions. The top most distributions in both panels correspond to the invariant 
yields of h measured at 40 degrees with the MRS spectrometer (scaled by 100 
for clarity purposes), followed by the yields of negative hadrons measured at 12 
(scaled up by 10) and 4 degrees respectively. More details about these distributions 
can be found in Ref. [ I13| 



System T^le l^Fits to power law shapes at n =r8.2. y 2 /NDF 

ar ) fit ar i meas 1 

GcV/c 

p + p 1.05 ± 0.06/0.95 ± 0.07 1.18 ± 0.16 10.9 ± 0.9 13. / 11 
d + Au 2.23 ± 0.09 / 2.1 ± 0.6 1.52 ± 0.1 12.3 ± 0.5 102. / 11 



2.2. Nuclear modification factor RdAu 

The d+Au system is compared to a "simpler" one: p + p where in this particular 
case, we do not expect the effects of saturation. This comparison is based on the 
assumption that the production of moderately high transverse momentum particles 
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scales with the number of binary collisions N co u in the initial stages. The so-called 
nuclear modification factor is defined as: 



RdAu 



1 N d Au(PT,V) 



Ncou N pp (p T ,ri) 

where N™u is estimated to be eaual to 7.2 ± 0.3 for minimum bias collisions. 

2 r 
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Fig. 2. Nuclear modification factor for charged hadrons at pseudorapiditics 
rj = 0,1.0,2.2,3.2. Statistical errors are shown with error bars. Systematic er- 
rors are shown with shaded boxes with widths set by the bin sizes. The shaded 
band around unity indicates the estimated error on the normalization to (N co u). 
Dashed lines at pr < 1 GeV/c show the normalized charged particle density ratio 

1 dN/dr){d+Au) 
(N coll ) dN/d V (pp) ■ 



Figure |21 shows the nuclear modification factor defined above for four 77 values. 
At mid-rapidity (77 = 0), the nuclear modification factor exceeds 1 for transverse 
momenta greater than 2 GeV/c in a similar way as the measurements performed 
by Cronin at lower energies 

A shift of one unit of rapidity is enough to make the Cronin type enhancement 
disappear, and further increases in 77 decrease even further the nuclear modification 
factor RdAu- Further details from this Figure [21 can be found in Ref. [I13|. 

We see a one to one correspondence between the RdAu values at low px and 
the ratio 7^ — r dN iif?i d ^ A ^ as demonstrated in Fig. |21 where that ratio is shown as 

{N co u) dN/dr){pp) ° 

dashed lines at pt < 1. 
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2.3. Centrality dependence R cp 

In the context of saturation, the suppression of the overall number of gluons in 
the Au wave function depends on a power of the number of participating nu- 
cleons Np^. t ; the higher this number, the stronger the suppression. We extract 
the number of participants from the multiplicity of the collision measured in the 
range | 77 |< 2.2. To study the dependence on centrality or number of participants 
three data samples of different centralitics were defined according to the multi- 
plicity of each event, and scaled histograms in transverse momentum were filled: 
N cen trai(PT) = jr — N _ 2 o%(pt) for events with multiplicities ranging from to 
20%. ^semi — central (p T ) = j^N 30 _ 50 %(p T ) for semi-central events with multi- 
plicities ranging from 30 to 50%, and finally, N periph (p T ) = j^-^N 60 _ SQ %(p T ) for 
peripheral events with multiplicities ranging from 60 to 80% with N co u values listed 
in table □ 

With these histograms, two ratios were constructed: R(?p tral = y""' r °'^| and 

rp semi— central Nsemi-central (pt) 

U CP - N vf „ ph {p T ) ■ 

Because these ratios are constructed with events from the same data run, many 
corrections cancel out. The only correction that was applied to these ratios is related 
to trigger inefficiencies that become important in peripheral events. The dominant 
systematic error in these ratios stems from the determination of the average number 
of binary collisions in each centrality data sample. This error is shown as a shaded 
band around 1 in Fig. 

Table 2. N part and N co u values extracted from HIJING calculations for d+Au 



Centrality 


N part (Au) 


N p art(d) 


Ncoll 


CentralO - 20% 


12.5 


1.96 


13.6 ± 0.3 


Semi — centralSO — 50% 


7.36 


1.79 


7.9 ± 0.4 


Peripheral®) - 80% 


3.16 


1.39 


3.3 ±0.4 



The four panels of figure show the central R^jp tral (filled symbols) and semi- 
central Ji s ^™ % - central (open symbols) ratios for the four 77 settings. At 77 = on the 
left-most panel, the central events yields are systematically higher than those of the 
semi-central events, and at the right-most panel (77 = 3.2) the trend is reversed; the 
yields of central events are ~ 60% lower than those from semi-central events at all 
values of transverse momenta. 



3. Discussion 



Figures [5] and |3] can be described in the following way: at 77 = the Cronin like en- 
hancement is produced by coherent multiple scatterings of deuteron valence quarks 
on a saturated Au (for x < 0.01). The height of the enhancement increases with 
the centrality of the collisions because the number of scatterings is higher. As one 
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Fig. 3. Central (full points) and semi-central (open points) R cp ratios (see text 
for details) at pseudorapidities 77 = 0,1.0,2.2,3.2. Systematic errors (~ 5%) are 
smaller than the symbols. The ratios at all pseudorapidities are calculated for the 
average charge - . 



moves to a higher rapidity y, the probability of gluons emission grows as P ~ a s y 
and additional corrections have to be included. These corrections can be written as: 
d £n±) = cts(2N — N 2 ) within saturated and non-linear systems or d ^i^ = cts^N 
in linear non-saturated systems. The variable N in these equations is related to 
the density of gluons in the nuclei. The linear term on the right side of these equa- 
tions describes the emission of gluons and the quadratic term represents interactions 
between gluons that reduce their numbers. The numerator of the nuclear modifica- 
tion factor RdAu is growing slowly with rapidity because of the taming effect of the 
quadratic term in the equation mentioned above, while the denominator continues 
to grow because the p+p system is more dilute. A similar suppression at all pt is 
seen in Figure [2] The strength of that suppression is proportional to the number of 
participant nucleons in the gold ion because it implies a higher number of gluons 
in the system making the effect of the quadratic term more and more important. 
This effect is also seen in Figure as the pseudorapidity of the detected particles 
changes from to 3.2 the suppression is stronger for the central sample of events. 

A recent calculation of the nuclear modification factor RdAu based on a de- 
scription of the Au wave function as a color glass condensate and the deuteron as 
a dilute system of valence quarks [ Ej agrees well with the BRAHMS results at 
i] = 3.2. 
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